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The  information  presented  in  this  report  was  generated  during  the 
performance  of  the  task  Oil  of  CONTRACT  No.  F33615-81-C-2050/ 
SB5448-81-C-0100.  The  technical  work  was  carried  out  at  the  University 
of  Missouri-Rolla,  Department  of  Mechanical  and  Aerospace  Engineering. 

Dr.  John  W.  Sheffield  was  the  principal  investigator  of  the 
program.  Raghav  P.  Kher  was  the  graduate  student  supported  by  this 
contract.  This  is  an  interim  report  which  presents  the  results 
generated  during  the  period  of  1  January  1984  to  1  October  1985.  The 
program  was  sponsored  by  the  Aeronautical  Systems  Division  of  the  Air 
Force  with  Ms.  V.  J.  Van  Griethuysen  of  the  Power  Technology  Branch  at 
the  Aero  Propulsion  Laboratory  (AFWAL/POOS)  serving  as  the  technical 
monitor. 
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I.  INTRODUCTION 

The  reliability  of  a  spacecraft  can  be  enhanced  greatly  if  its  dependence 
upon  power  and  endurance  of  moving  parts  is  minimized.  In  reviewing  the 
projections  of  future  spacecraft  power,  life  and  mission  requirements  that  are 
currently  established,  the  thermal  management  of  spacecraft  is  seen  as  an 
important  issue  in  the  overall  design  of  the  spacecraft  system.  Future 
spacecraft  thermal  management  systems  (SCTMS)  will  require  advanced 
technologies  in  order  to  achieve  the  desired  weight,  size  and  life  goals. 

These  technologies  might  include  high-capacity  heat  pipes,  phase  change 
materials  (PCMs',  and  deployable  radiator  panels.  High-capacity  heat  pipe 
technologies  offer  opportunities  to  customize  isothermal  panel  design  and 
allow  deployment  of  additional  thermal  control  surfaces.  The  emerging 
two-phase  heat  pipe  technologies  [1]  include  axial  groove,  monogroove, 
capillary  pumped,  and  mechanically  pumped  loop  configurations.  The  heat 
transport  capacity  of  these  heat  pipes  are  shown  in  Table  I. 

Table  I.  Emerging  Heat  Pipe  Technologies 
TECHNOLOGY  CAPACITY  [kW  m] 

Axial  Groove  (current)  0.25  -  1.0 

Monogroove  ( current/projected)  25.0  -  100.0 

Capillary  Pumped  (projected)  250.0  -  1000 

Mechanically  Pumped  (projected)  25,000.0 

Reviewing  earlier  developments  of  spacecraft  technology  such  as  the 
efforts  sponsored  by  NASA  (National  Aeronautics  and  Space  Administration)  to 
develop  PCM  thermal  control  devices,  passive  thermal  control  systems  were 
shown  to  be  feasible  and  test  results  showed  that  PCM  systems  could  reduce 
onboard  equipment  temperature  fluctuations  by  at  least  75%  compared  to 
unprotected  spacecraft  systems  [2].  Basically,  this  earlier  passive  thermal 
management  system  concept  was  to  sandwich  the  PCM  between  the  equipment  to  be 


> _ 


£%*> 
i  . 


i 


controlled  and  the  outer  radiating  surface  of  the  satellite.  A  recent  paper 
[3]  investigated  the  sensitivity  of  the  future  SCTMS  weight  to  the  choice  of 
PCM.  The  candidate  PCMs  can  be  classified  as  form  stable  polymers, 
encapsulated  salt  hydrates,  paraffins  and  eutectic  alloys. 

Thermal  management  of  future  spacecraft  will  introduce  many  challenging 
design  requirements.  The  thermal  control,  heat  storage,  heat  transport,  and 
heat  rejection  represent  a  set  of  major  design  problems  when  the  mission 
constraints  of  large  thermal  loads  are  included.  One  of  the  more  challenging 
thermal  design  requirements  is  the  capability  to  handle  both  large  and  pulsed 
thermal  loads.  During  and  immediately  following  burst  power  conditions,  large 
thermal  loads  must  be  removed  from  the  source  and  rejected  in  some  manner  that 
minimizes  the  overall  burden  on  the  total  system.  The  application  of  PCM 
appears  to  be  an  attractive  choice  [4].  The  design  of  a  SCTMS  requires  the 
determination  of  a  large  number  of  variables  and  their  interdependence, 
including  the  integration  of  both  mission  and  launch  payload  requirements. 

The  comparison  of  conceptual  designs  requires  system  studies  of  configuration, 
thermal  structures  and  subsystems.  Figure  1  gives  an  example  of  the  philosophy 
of  design  concepts  for  SCTMS.  Figure  2  shows  the  steps  in  the  overall  thermal 
control  design  process  [5]. 

It  has  become  clear  to  the  designers  of  SCTMS  that  significant  reduction 
in  the  overall  weight  of  space  power  systems  can  be  achieved  through  the  use 
of  advanced  radiator  concepts,  such  as  modular  heat  pipe  panels,  pumped  fluid 
loop,  deployable  and  liquid  droplet.  The  variable  conductance  heat  pipe 
thermal  control  systems  can  provide  near  isothermal  conditions  without  moving 
parts  and  are  finding  increasing  acceptance.  The  critical  element  in  a 
deployable  radiator  design  is  the  thermal  contact  of  the  joints.  A  flexible 
heat  pipe  is  an  excellent  candidate  for  this  application. 
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Figure  1.  Philosophy  of  Design  Concepts 


Recently,  as  a  part  of  an  investigation  'Weight  Characteristics  of  Future 
Spacecraft  Thermal  Management  Systems'  [6]  sponsored  by  Air  Force  Wright 
Aeronautical  Laboratories'  Power  Technology  Branch,  a  SCTMS  concept 
incorporating  a  PCM  into  a  variable  conductance  heat  pipe  was  explored.  This 
concept  of  PCM  microspheres  allows  for  maximized  surface  to  volume  ratio  and 
very  rapid  charging  and  discharging  of  thermal  load.  The  microspheres  have  an 
additional  potential  of  acting  as  a  wicking  material  for  condenser  of  heat 
pipes  making  this  concept  very  promising. 

The  weight  characteristics  of  future  spacecraft  thermal  management 
systems  have  been  discussed  in  a  recent  paper  [3].  Performance  parameters 
include  the  peak  thermal  load,  the  peak  to  average  thermal  load  ratio  and  the 
length  of  the  peak  thermal  load.  In  addition,  the  sensitivity  of  SCTMS  weight 
to  the  choice  of  PCM  was  determined  via  a  heat  of  fusion  parametric  analysis. 
As  an  ongoing  research  task,  this  report  includes  extensive  literature  survey 
covering  different  concepts  for  SCTMS,  and  hardware  requirements  of  their 
applications.  The  first  few  pages  of  this  report  summarizes  different  aspects 
considered  before  choosing  the  heat  pipe/PCM  concept  for  current 
investigation.  Detailed  laboratory  tests  on  different  PCMs  were  done  and  a 
latent  thermal  energy  storage  unit  was  built  and  tests  were  carried  out  to 
determine  the  performance  of  the  unit. 

To  assist  in  the  design  process,  a  knowledge-based  system  has  been 
developed  in  a  current  artificial  intelligence  programming  language,  PROLOG 
(Programming  in  Logic)  [7].  This  software  serves  as  a  design  assistant  for 
the  development  of  SCTMS.  It  supports  facilities  such  as  database  retrieval, 
numeric  computation  and  logical  decision  making.  These  features  allow  one  to 
automatically  consider  all  relevant  combinations  of  components,  technologies 
and  materials  for  the  design  problem. 


II.  BACKGROUND 


A.  PHASE  CHANGE  MATERIALS 

A  promising  thermal  control  concept  not  yet  widely  employed  in  the  space 
is  the  use  of  Phase  Change  Materials  (PCMs).  These  materials  have  extremely 
long  lifetime,  and  could  be  used  to  maintain  a  relatively  constant  temperature 
of  an  electronic  device  or  to  store  heat  energy  for  later  release.  The  latent 
heat-of-fusion  concept,  that  is  storing  and  recovering  heat  through  the 
solid-liquid  phase  change  process,  has  two  important  advantages:  (1)  the 
latent  heat  of  fusion  of  most  materials  is  much  higher  than  their  sensible 
heat,  thus  requiring  much  smaller  mass  of  storage  medium  for 
storing/recovering  a  given  quantity  of  thermal  energy  and  (2)  the  thermal 
storage  process  occurs  at  nearly  constant  temperature,  which  is  typically 
desirable  for  efficient  operation  of  most  thermal  systems.  A  good 
understanding  of  the  heat  transfer  process  involved  is  essential  for 
accurately  predicting  the  thermal  performance  of  the  system  and  for  avoiding 
costly  system  design. 

The  heat  transfer  and  design  aspects  of  Latent  Heat  Thermal  Energy 
Storage  (LHTES)  systems  are  more  complex  than  those  in  sensible  systems  due 
to:  (1)  nonlinearity  of  the  problem  resulting  from  the  motion  of  the 
solid-liquid  interface  during  phase  change  process;  (2)  uncertainity  of  the 
interface  thermal  resistance  between  the  container  wall  (or  transport  fluid 
channel)  and  the  solid  PCM;  (3)  volume  change  of  phase  change  (upon  shrinkage, 
void  cavities  may  arise  in  the  solid  and  reduce  the  heat  transfer  process 
significantly);  and  (4)  the  presence  and  configuration  of  voids  in  solids 
which  may  also  introduce  the  complex  phenomena  of  natural  convection  and 
internal  radiation,  particularly  at  higher  temperatures .  There  are  few 
examples  of  technical  data  available  on  the  latter  four  problem  areas 
mentioned  that  would  be  applicable  to  the  design  of  latent  heat  thermal  energy 
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storage  systems  (LHTES).  The  following  pages  give  a  few  design  concepts  for 
LHTES.  Table  II  gives  the  desired  characteristics  of  latent  heat  thermal 
storage  systems. 

Table  II.  Desired  Characteristics  of  Latent  Thermal  Storage  System 

1.  Compact:  Large  heat  storage  capacity  per  unit  mass  volume. 

2.  High  storage  efficiency:  High  thermal  conductivity. 

3.  Heat  storage  medium  with  suitable  properties  in  operating 
temperature  range. 

4.  Capacity  to  charge  and  discharge  with  the  largest  heat 
input/output  rates  but  without  large  temperature  gradients. 

5.  Ability  to  undergo  large  number  of  charging/discharging  cycles 
without  loss  in  performance  and  storage  capacity. 

6.  Small  self-discharging  rate  i.e.,  negligible  heat  losses  to 
surroundings. 

7.  Inexpensive 

1.  Micro-encapsulated  PCM  Concept:  In  this  concept,  relatively  large 
quantities  of  PCM  are  encapsulated  with  thin  layers  of  containment  material 
(steel)  which  deform  to  accommodate  volumetric  changes  associated  with  change 
of  phase  [8].  Relatively  small  quantities  of  containment  material  would  be 
required,  since  the  thickness  would  have  to  be  small  to  minimize  bending 
stress.  One  such  configuration,  consisting  of  long  'plank1  shaped  blocks  of 
PCM  ib  shown  in  Fig.  3.  The  plank  cross  section  would  be  essentially 
rectangular  when  the  PCM  is  at  maximum  density  in  solid  state.  The  edges 
would  be  relatively  stiff  and  the  top  and  bottom  surfaces  would  bulge  outward 
as  indicated  to  accommodate  volumetric  expansion  associated  with  melting.  The 
planks  would  be  stacked  in  a  pressure  vessel  as  shown  in  Fig.  3-B,  with  the 
gap  regions  provided  for  liquid  heat  transfer  fluid  flow.  The  gap  height 


should  be  adjusted  by  design  of  the  plank  edge  shape  to  be  as  small  as 
possible  consistent  with  allowable  flow  pressure  drops  and  required  film 
coefficient  for  heat  transfer.  High  overall  packing  density  is  desirable  to 
minimize  pressure  vessel  size  and  cost. 

The  principal  disadvantage  of  this  concept  is  the  cost  of  the  required 
pressure  vessels.  An  intermediate  low  pressure  circulating  loop  might  be 
considered  to  eliminate  the  pressure  vessels,  but  the  cost  of  the  heat 
exchangers  for  such  a  loop  would  probably  be  comparable  to  the  cost  of  the 
equivalent  boiler  capacity  being  replaced.  A  better  approach  primarily  suited 
to  very  large  thermal  energy  storage  applications. 

2.  Intermediate  Pumped  Loop:  The  use  of  a  low  pressure  intermediate 
liquid  metal  circulating  is  an  alternative  variation  on  the  base  line 
tube/shell  PCM  concept.  The  low  pressure  would  permit  significant  reductions 
in  wall  thickness  for  the  tubes  and  shell  domes.  This  would  be  off-set  to  some 
extent  by  the  cost  of  heat  exchangers,  pumps,  and  piping  for  the  liquid  metal 
loop.  Heat  exchanger  bypass  lines  could  be  provided  in  the  loop  with  bypass 
flow  varied  to  compensate  for  variable  thermal  resistance  within  the  PCM 
associated  with  movement  of  liquid  solid  interface.  This  bypass  control 
should  result  in  more  uniform  steam  and  water  outlet  temperatures,  however,  it 
might  experience  reduce  reliability. 

3.  Fluidized  Encapsulated  PCM:  Micro-encapsulating  the  PCM,  using 
particle  sizes  comparable  to  that  of  pulverized  coal  (e.g.  50 - 1 0 0 p)  is 
attractive  because  it  eliminates  the  need  for  heat  exchangers  though  this  may 
be  off-set  by  encapsulation  cost  [8].  This  concept  is  illustrated  in  Fig.  4 
using  a  one  bin  empty  system.  For  heat  storage,  cold  encapsulated  PCM  is 
drawn  from  a  filled  bin,  entrained  by  carrier  steam,  and  deposited  hot  in  the 
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figure  3.  Microencapsulated  PCM  a)  Cross-Section 

b)  Stacking  Arrangement 


next  (empty)  bin.  One  bin  is  depleted  as  the  other  fills.  When  full,  valves 
are  adjusted  to  repeat  the  process  of  taking  PCM  from  another  cold  bin  and 
putting  it  into  the  previously  emptied  bin.  Bins  would  be  large  enough  to 
permit  particle  settling,  with  the  carrier  fluid  drawn  off  at  the  top.  For 
heat  removal,  the  carrier  fluid  would  be  liquid  water  in  the  feedwater  heating 
case,  or  liquid  evolving  to  steam  in  auxiliary  power  cycle  case.  As  with  the 
other  PCM  concepts,  several  PCMs  with  different  melt  temperatures  might  be 
used,  using  separate  sets  of  bins  for  each  PCM,  each  based  on  one  bin  empty 
concept.  The  micro-encapsulated  PCM  -  system  would  differ  from  a  pulverized 
coal  system  in  having  to  take  full  cycle  pressure.  Provision  for  bin-to-bin 
venting  would  be  required  to  minimize  pressure  variations  during  PCM  transfer. 
Also,  it  would  be  desirable  to  have  following  mass  flow  rate: 

mw  (hsw)  mp  (  H  +  hsp) 
where  hsw  =  sensible  heat  of  water 

hsp  =  sensible  heat  of  PCM 

mw  =  mass  of  water 

so  that  all  the  heat  transfer  for  each  PCM  could  be  accomplished  with  a  single 
entrainment.  Then  sensible  energy  storage  would  be  greater  for 
microencapsulated  PCM  than  for  other  concepts,  since  the  PCM  would  quickly 
approach  the  carrier  fluid  temperature  with  negligible  exit  temperature 
differences.  This  permits  smaller  quantities  of  PCM  while  all  systems  will 
have  containment  material  associated  with  sensible  heat  storage,  the  mass  of 
such  material  may  be  somewhat  larger  for  microencapsulated  than  it  would  be 
for  other  concepts. 

From  a  safety  point  of  view,  microencapsulation  presents  two  problems, 
failure  of  encapsulating  shells  and  failure  of  the  pressurized  bins. 
Encapsulated  shell  failure  would  release  small  quantities  of  PCM  which  could 
react  with  carrier  fluid,  with  some  particle  aglomeration  probable.  Use  of 
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relatively  course  mesh  screen  or  periodic  bin  cleaning  might  be  required.  Bin 
structural  failure,  as  with  any  pressure  vessel,  is  a  serious  problem,  but 
only  the  void  volume  in  the  bin  would  contain  high  pressure  fluid  capable  of 
contributing  work  of  expansion.  Also,  the  encapsulated  particles  should 
present  somewhat  less  hazard  than  bulk  PCM. 

Major  drawbacks  associated  with  materials  and  processes  for  encapsulation 
of  microencapsulated  are  not  completely  known.  If  encapsulated  PCM  could  be 
provided  at  moderate  cost  ($0 .5/Kg) ,  the  system  would  be  more  attractive  but 
would  still  result  in  added  cost  for  pressure  vessel  and  high  pressure  valves 
[8].  Development  of  hardware  and  flow  relations  associated  with  handling  of 
microencapsulated  PCM  are  not  trivial  problems  but  appear  soluble  with 
straight  forward  engineering  effort. 

4.  Moving  PCM  Concept:  The  moving  PCM  concept  attempts  to  minimize 
required  heat  transfer  surface  area  by  eliminating  most  of  the  thermal 
resistance  within  the  PCM  [8].  The  liquid  solid  interface  is  always  kept 
close  to  the  heat  transfer  surface.  This  concept  is  illustrated  schematically 
in  Fig.  5.  During  usage,  liquid  PCM  flow  through  slit  nozzles  on  to  rotating 
drums.  The  heat  sink  fluid  flows  through  the  drums,  cooling  the  surfaces  and 
freezing  the  PCM.  The  solid  PCM  is  scrapped  off  after  making  a  partial 
rotation,  falling  into  a  storage  bin.  For  heat  storage,  the  source  fluid  is 
fed  through  heat  exchangers  located  at  the  bottom  of  the  bin  which  has  a 
pitched  floor  arranged  so  that  the  PCM  melts  off  as  it  drains.  The  liquid  PCM 
is  pumped  into  an  insulated  tank  until  needed  for  part  of  the  cycle.  Potential 
problems  associated  with  liquid  PCM  pumping  are  recognized,  such  as  flow 
blockage  in  the  pumps  due  to  PCM  solidification.  This  concept  was  originally 
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conceived  as  a  terrestrial  design,  however,  a  closed  loop  pumped  system  could 
be  considered  for  space  designs. 

5.  Heat  Pipe  Concept  :  This  is  a  new  type  of  a  passive  heat  exchanger 
presently  under  development.  The  major  advantages  offered  by  this  type  of 
heat  exchanger  are  its  modular  design  and  flexibility  in  application  in  that 
it  can  equally  be  used  with  liquid  or  air  as  heat  transfer  medium. 

Figure  6  shows  a  schematic  of  a  single  heat  exchanger  module  essentially 
compromised  of  a  container  (1)  of  square  cross  section  provided  with  a  heat 
pipe  (2)  along  its  longitudinal  axis  [9].  The  container  as  well  as  heat  pipe 
are  divided  into  three  regions  -  region  A,  B  and  C  in  Fig.  6  by  separation 
walls  (5).  One  of  the  regions  B,  is  filled  with  a  heat  of  fusion  type 
material  and  is  called  storage  chamber.  The  remaining  two  regions  A  and  C  are 
in  contact  with  fluid  flowing  through  the  heating  loops,  and  are  termed  as 
heat  source  or  heat  sink  regions  respectively. 

The  heat  pipe  length  within  the  storage  chamber  is  provided  with  closely 
and  equally  spaced  square  fins  of  large  cross  sections  (3),  made  from  a 
material  of  high  thermal  conductivity,  i.e.  aluminum.  The  storage  medium  (4) 
fills  the  free  volume  between  the  fins.  Thus  embedded  within  the  storage 
medium,  the  fins  compensate  for  its  generally  poor  thermal  conductivity  and 
simultaneously  insure  small  heat  flow  paths  within  LHTES  system. 

The  use  of  the  heat  pipe  offers  several  advantages  and  increases 
flexibility  in  operation  and  application.  Some  of  these  are: 

a.  Additional  heat  exchangers  are  eliminated.  Furthermore,  the  heat  pipe 
length  within  the  heat  source  and  heat  sink  regions  may  be  suitably  finned 
depending  on  whether  air  or  liquid  is  employed  as  heat  transfer  medium  in 
these  regions. 


b.  The  heat  pipe  transports  heat  under  low  temperature  gradient  so  that 
almost  isothermal  heat  source  in  contact  with  the  fins  and  thermal  storage 
medium  within  the  storage  chamber  is  attained. 

c.  The  heat  flux  transformation  capability  of  the  heat  pipe  can  be 
utilized  to  give  low  heat  flux  densities  within  the  storage  chamber  for  large 
heat  flow  rates  in  heat  source/sink  sections. 

d.  The  heat  pipe  can  operate  unidirectionally  as  a  thermal  diode. 

B.  HEAT  PIPES 

The  heat  pipe  is  a  device  with  an  effective  thermal  conductance  much 
higher  than  even  the  best  heat  conducting  metal.  This  property  makes  it 
especially  suited  for  space  applications  in  which  a  large  amount  of  heat  must 
be  transferred  from  one  area  to  another,  or  when  a  stable  'isothermal' 
condition  is  required  in  an  area  subject  to  internal  heat  generation. 

Figure  7  shows  the  basic  concept  for  a  standard  heat  pipe  [5],  The  heat 
pipe  operates  on  the  same  principle  as  a  reflux  boiler,  but  depends  upon  the 
capillary  action  of  a  wick  rather  than  gravity  to  return  the  condensed  vapor 
from  the  heat  radiator  (condenser)  to  the  heat  source  (evaporator). 

The  standard  heat  pipe  has  a  fixed  conductance  and  must  therefore  be 
designed  to  handle  a  constant  and  specific  heat  load.  A  later  development, 
the  variable  conductance  heat  pipe,  makes  it  possible  to  control  the 
conductance  as  a  function  of  temperature,  which  broadens  considerably  the 
applications  for  this  highly  efficient  device.  Variable  conductance  is 
achieved  by  means  of  a  noncondensable  gas  contained  at  the  condenser  end  of 
the  heat  pipe.  This  gas  expands  and  contracts  with  changes  in  vapor 
temperature  within  the  pipe,  thus  changing  the  effective  condenser  heat 
transfer  area. 


C.  RADIATORS 


In  almost  all  advanced  space  projects,  there  are  requirements  to  remove 
heat  from  space  vehicles  during  very  long  periods  of  time.  The  only  practical 
way  of  such  continuous  heat  rejection  from  a  system  in  an  outer  space 
environment  is  by  radiation. 

Radiators  for  space  vehicle  environmental  control  systems  may  be 
considered  as  low-temperature  radiators;  they  are  differentiated  from 
high-temperature  radiators  in  that  the  net  heat  rejection  from  low-temperature 
radiators  is  of  the  same  order  of  magnitude  as  the  heat  flux  incident  upon  the 
radiator  from  surrounding  environment.  In  other  words  the  radiator 
temperature  is  relatively  low  and  is  not  greatly  above  the  effective 
environment  temperature.  The  distinguishing  feature  of  low  temperature 
radiators  is  that  they  depend  to  a  large  extent  upon  the  spectrally-selective 
material  used  as  a  surface  coating.  This  coating  should  have  as  high  as  an 
emissivity  as  possible  and  as  low  as  possible  absorptivity  for  solar 
radiation.  The  design  of  radiators  requires  knowledge  of  the  heat  transfer 
processes  and  thermal  design,  methods  of  meteoroid  protection,  fabrication, 
and  reliability  requirements  for  radiator  components.  Although  the  basic  laws 
of  radiant  heat  transfer  are  well  known,  the  design  of  the  actual  hardware  to 
be  used  in  radiator  systems  for  space  vehicles  is  not  well  developed.  One 
reason  for  this  is  that  there  are  many  severe  criteria  that  must  be  met  by 
space  heat  rejection  systems.  Among  these  criteria  are  low  weight,  high 
reliability,  operation  in  vacuum  environment  of  space,  and  the  capability  of 
continuous  operation  despite  many  possible  collisions  with  meteoroid.  The  low 
weight  and  invulnerability  to  meteoroid  collision  are  probably  the  most 
limiting  specifications  to  be  considered  in  radiator  design. 


0.  ARTIFICIAL  INTELLIGENCE 


Definitions  of  the  phrase  'artificial  intelligence*  (AI)  vary  widely  from 
author  to  author,  and  a  given  author  will  tailor  his  or  her  definitions  to  his 
or  her  particular  audience.  Definitions  made  by  AI  researchers  to  their 
colleagues  emphasize  subtle  points  of  philosophy,  technology,  and  methodology. 
Definitions  for  a  more  general  audience  must  give  a  more  general 
characterization.  Many  definitions  refer  to  the  long  term  goals  of  artificial 
intelligence  research  to  emulate  human  intellectual  abilities.  This  may  give 
a  philosophical  perspective  on  the  domain  but  it  may  leave  the  impression  that 
accounting  and  inventory  programs  are  really  artificial  intelligence  since 
they  emulate,  and  in  fact  replace,  humans  that  used  to  perform  these 
functions.  Hence,  it  is  useful  to  identify  the  types  of  problems  which  have 
been  the  subject  of  AI  research. 

A  recent  paper  [10]  identifies  three  features  of  typical  AI  problems. 
These  problems  have  been  characterized,  first  of  all,  as  combinatorial ly 
complex  problems.  The  second  feature  of  AI  problems  is  that  AI  has  been 
applied  to  problems  for  which  no  known  'solution  process'  exists.  The  only 
path  to  a  solution  may  be  through  the  application  of  'general  principles', 
'rules  of  thumb',  or  'educated  guesses'.  The  problems  are  computationally  so 
complex  that  it  is  impossible  to  ever  find  an  optimal  solution.  In  fact,  the 
problems  involve  the  consideration  of  so  many  factors  that  it  may  be 
impossible  numerically  to  quantify  the  'value'  of  a  solution.  It  may  only  be 
possible  to  qualify  a  solution  as  good  or  bad,  or  qualify  one  solution  as 
better  or  worse  than  another.  Typically,  AI  problems  are  so  hard,  that  just 
one  satisfactory  solution  is  valued. 

AI  is  distinguished  by  the  problems  that  it  addresses  and  also  the 
techniques  and  technology  which  it  applies  to  those  problems.  Kempf  [11] 
defines  AI  as  a  three  stage  problem  solving  technique  to  determine: 


a.  All  available  knowledge  that  can  be  applied  to  the  combinatorial 

complexity  of  the  problem. 

b.  Some  representation  for  that  knowledge. 

c.  A  methodology  for  applying  that  knowledge  to  the  problem. 

Of  course  this  is  a  good  problem  solving  methodology  for  any  subject  domain. 
For  familiar,  well  defined  engineering  problems,  the  relevant  knowledge  is 
applied  to  equations  and  algorithms  where  the  representation  is  often  a 
FORTRAN  program.  This  knowledge  is  then  utilized  by  applying  the  program  to  a 
file  of  input  data  yielding  an  output  file  of  solutions.  However,  as  problems 
become  more  complex  and  less  well  defined,  the  relevant  knowledge  covers  a 
spectrum  of  facts,  decision  rules,  equations,  and  algorithms.  As  problems  of 
representation  become  more  complex,  the  problem  of  applying  the  knowledge 
becomes  more  difficult.  For  this  reason  most  AI  applications  have  been 
implemented  in  symbolic  languages,  such  as  LISP  or  PROLOG,  which  are  more 
suitable  for  the  representation  of  many  different  kinds  of  information  and 
which  can  access  and  process  that  information  in  many  different  ways. 

The  first  goal  is  to  capture  the  expertise  of  a  human  researcher  [12], 
and  to  emulate  a  methodology  in  the  design  of  future  SCTMS.  When  all  of  the 
potential  design  options  are  considered  together  in  all  of  their  combinations 
the  number  of  potential  designs  becomes  exponentially  large.  Instead  the 
design  process  is  far  from  that  which  we  would  call  algorithmic.  It  involves 
a  creative  process  of  inventing  a  potential  design,  which  is  then  further 
refined  and  modified  according  to  many  general  principles.  Finally,  the 
solution  must  meet  certain  basic  requirements.  It  may  even  be  possible  to 
characterize  a  design  as  good  or  bad,  but  so  many  factors  must  be  considered 
that  it  is  not  clear  what  'optimal'  might  mean. 

Attempts  to  automate  design  processes  is  a  relatively  new  endeavor  in  AI . 
Many  other  areas  of  AI  have  been  the  subject  of  years  of  ongoing  research: 


game  playing,  theorem  proving,  vision,  natural  language  understanding, 
learning,  diagnosis,  etc.,  but  engineering  design  poses  some  unique 
difficulties.  A  human  designer  brings  a  vast  body  of  factual  and  experiential 
knowledge  to  bear  on  a  design  problem.  The  resulting  design  must  of  course 
satisfy  certain  requirements,  but  the  number  of  satisfactory  designs  is 
unlimited.  Hence,  an  AI  system  which  emulates  a  human  designer  must  have 
access  to  vast  bodies  of  knowledge,  it  must  have  the  capability  to  learn  from 
its  previous  designs,  and  it  must  have  very  sophisticated  mechanisms  for 
focusing  on  a  few  candidate  designs. 

A  recent  work  on  automated  design  [10]  gives  some  direction  to 
implementation  of  automated  design  systems.  The  authors  of  that  paper 
identify  three  classes  of  design  problems: 

a.  Class  1  is  very  creative,  requires  extensive  knowledge  and 
experience,  and  involves  consideration  of  an  unlimited  number  of  solutions. 

b.  Class  2  is  generally  routine,  perhaps  even  repetitive,  but  it 
requires  in  most  cases  some  inventiveness  and  creativity  to  resolve  some 
aspects  of  the  design. 

c.  Class  3  follows  a  well-established,  reasonably  well  understood  design 
procedure.  However  the  design  problem  is  still  'hard'  in  the  sense  that  there 
is  still  an  effectively  unlimited  number  of  potential  designs.  Very  specific 
knowledge  and  expertise  is  required  to  focus  on  a  few  candidate  designs  and 
develop  them  from  initial  concept  through  several  stages  of  modification  and 
refine  ",r 

The  authors  of  that  paper  proposed  that  the  Class  3  design  is  within  the 
reach  of  current  AI  technology  and  techniques.  Since  the  design  of  future 
SCTMS  is  hardly  a  routine  process,  the  first  step  of  this  research  was  to 
restrict  the  scope  of  the  design  problem.  The  work  effort  was  focussed  on  a 
specific  configuration  (heat  pipe,  heat  storage,  radiator)  where  a  system 


which  would  substantially  automate  the  process  of  selecting  compatible 
components  for  the  system,  properly  sizing  them,  and  refining  their 
characteristics  was  implemented. 


III.  PROPOSED  CONCEPT 

In  this  investigation  a  simplistic  approach  was  used  in  the  modeling  of 
the  SCTMS.  The  basic  concept  was  to  divide  the  system  into  the  generic 
components  for  the  thermal  modelling.  Figure  8  shows  some  design  aspects  to 
the  proposed  concept  formulation.  Figures  9a  and  9b  illustrated  schematically 
the  baseline  concept.  Three  components  were  selected.  They  are  the  heat 
transport  heat  pipes,  the  thermal  storage  system  using  a  phase  change  material 
and  the  space  radiators.  A  transient  thermal  analysis  was  conducted  to 
determine  the  response  for  pulsed  thermal  loads.  The  transient  thermal 
analysis  was  the  Lumped-System  Analysis  [13]  where  the  spatial  variation  of 
temperature  is  neglected  and  the  variation  of  the  temperature  of  the  component 
with  respect  to  time  is  studied.  In  this  method  the  geometry  of  the  component 
is  immaterial  since  the  temperature  is  considered  to  be  a  function  of  time 
I  only,  hence  the  analysis  becomes  very  simple.  A  three-lump  system  was  used  to 

extend  the  study  of  the  transients  in  a  composite  system  consisting  of  the 
three  different  compoents.  In  the  following  analysis,  we  considered  the 
j  three  components  of  the  SCTMS,  i.e.,  heat  pipe,  PCM  and  radiator,  to  be  lump 

masses  having  essentially  uniform  temperatures. 

A.  MODELING  OF  THE  PROPOSED  SYSTEM 


1.  Transient  Thermal  Analysis:  The  mathematical  formulation  of  the 
problem  depends  on  the  number  of  lumps  considered  for  the  system  analysis  and 
on  the  type  of  boundary  conditions  for  the  heat  transfer.  In  this  analysis, 


Figure  9.  Baseline  Concept  a)  Conceptual  Heat  Rejection  System 


heat  transfer  boundary  conditions  of  the  third  kind  were  used  at  all  of  the 
boundaries.  The  energy  equation  for  a  component  may  be  stated  as  the 
following: 


The  net  rate  of  heat  flow  into  a  component  through  the  boundaries  is 
equal  to  the  rate  of  increase  of  the  internal  energy  of  the  component. 

For  the  three-lump  system,  the  energy  balances  for  the  three  components 
respectively  are  as  follows: 


h  dVdt 

-  Q  “  Qh 

(HEAT  PIPE) 

[1] 

Hp  dX/dt 

^h  “  Qp 

(PCM) 

[2] 

Cr  dTr/dt 

^p  "  ^r 

(RADIATOR) 

[3] 

Here  M  is  the  mass  of  the  component,  C  is  the  specific  heat  of  the  component, 

T  is  the  temperature  of  the  component,  X  is  the  melt  fraction  of  the  phase 
change  material,  t  is  the  time  variable,  H  is  the  latent  heat  of  fusion  of  the 
phase  change  material,  and  the  Q's  are  the  heat  flow  rates  of  the  three 
components.  The  three  heat  flow  rates  are  defined  as  follows: 


yt)  =  Uh  Ah  (  Th(t)  -  Tp)  [4] 

Qr(t)  =  Up  Ap  (  Tp  -  Tr(t)  )  [5] 

Qp(t)  =  Ur  Ar  (  Tr(t)  -  T$  }  [6] 

Here  the  ' U *  parameters  are  the  overall  heat  transfer  coefficients  for  the 
three  components,  and  the  'A'  parameters  are  the  associated  areas  of  heat 
transfer  of  the  three  system  components.  Substituting  Eqs.  [4-6]  into  Eqs. 
[1-3]  results  in  three  linear,  first  order,  ordinary  differential  equations 
for  the  transient  variations  of  the  temperature  of  the  heat  pipes,  the  melt 
fraction  of  the  phase  change  material  and  the  temperature  of  the  radiators. 

The  heat  flow  from  the  thermal  bus  is  assumed  to  be  a  periodic  function  of 
time.  For  simplicity,  a  square  wave  function  was  selected  as  the  periodic 
function.  It  can  be  shown  that  for  the  heat  flow  represented  by  a  square  wave 
that  the  ratio  of  the  'peak-to-average '  heat  flows  is  equal  to  one  plus  the 
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ratio  of  the  'off-to-on'  time  periods.  For  example,  the  'peak-to-average' 
heat  flow  ratio  for  the  special  case  of  equal  'on'  and  'off'  is  equal  to  two. 

The  reference  masses  of  the  three  components  of  the  thermal  management 
system  for  the  spacecraft  are  defined  by  the  following  expressions: 


Mh  -  Qh  / 

Gh 

[7] 

Mp  =  Qp  n 

tp  /  H 

[8] 

Mr  =  Qr  / 

Gr 

[9] 

Here  the  parameter  ' tp ’  is  the  time  of  the  peak  heat  load,  the  parameter  'n' 
is  the  number  of  peak  time  intervals,  and  the  ' G*  parameters  are  the  specific 
weights  of  the  components. 

2.  Thermal  Management  Subsystems'  Models:  The  specific  weights  of  the 
heat  pipes  and  the  radiators  are  evaluated  at  a  peak-to-average  heat  flow  of 
one,  i.e.  constant  heat  flow  for  the  reference  values.  The  effective, 
time-averaged,  specific  weights  of  the  heat  pipes  and  radiators  can  be 
calculated  once  the  actual  heat  flows  are  obtained  from  the  solutions  of  Eqs. 
[1-3].  The  time-averaged  specific  weight  of  a  heat  pipe  can  be  defined  as  the 
time-averaged  heat  flow  divided  by  the  mass  of  the  heat  pipe,  and  hence  can  be 
defined  by  the  following  expression: 

Sh  *  °Ah'"h  ■  (  <V  Qft)  sh  no] 

Similarly,  the  time-averaged  specific  weight  of  a  radiator  can  be  defined  as 
the  time-averaged  heat  flow  divided  by  the  mass  of  the  radiator,  hence  can  be 
defined  by  the  following  expression: 

G#r=  0r/  »r-  (  0Ar/  0A)  Gr  [11] 

The  time-averaging  is  simply  the  time  integral  of  the  function  divided  by  the 
time  interval.  The  total  mass  of  the  thermal  management  system  is  the  sum  of 
the  three  components:  the  heat  pipes,  the  PCM  and  the  radiators. 


Masses  of  the  SCTMS  Components 


a.  Heat  Pipe;  To  determine  the  masses  of  the  heat  pipes  in  a  SCTMS,  a 
specific  mass  formulation  is  used. 

d"h  -  '  «h  [12 

where  Nh  is  the  number  of  heat  pipes. 

The  total  number  of  heat  pipes  can  be  calculated  from  the  total  heat  load 
divided  by  the  heat  load  of  a  single  heat  pipe,  i.e. 

dNh  ■  dOt  /  Qh  [13 

The  heat  load  of  a  single  heat  pipe  can  be  calculated  from  the  maximum  heat 
pipe  capacity  divided  by  the  maximum  heat  pipe  length,  i.e. 

°h  ■  CAh  1  Lh  [14 

where  CA^  =  Capacity  of  the  heat  pipe. 

Lh  =  Length  of  the  heat  pipe. 

Combining  Eqs.  12  and  13  yields 

dMh  =  dot  /  ch  [is; 

The  evaporator  heat  fluxes  have  constraints  given  by 

qe  -  Qh  /  (Le  »  oh)  [ie; 

Thus  combining  with  Eq.  14  yields  the  evaporator  length 

Le  ’  CAh  '  lie  Lh  w  °h)  [11 

where  Lg  =  Length  of  the  evaporator. 

Dg  =  Diameter  of  the  heat  pipe. 
qe  =  Heat  flux  of  evaporator. 

Similarly  the  condenser  heat  fluxes  have  limitations  given  by 

qc  -  Qh  /  (Lc  71  °h)  CIS 

Thus  combining  with  Eq.  14  yields  the  condenser  length 

Lc  -  CA„  /  (qc  Lh  »  0h)  [19 

qc  =  Heat  flux  at  the  condenser. 

1.  =  Length  of  the  condenser. 
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b.  Phase  Change  Material:  To  determine  the  weight  of  phase  change 
material  in  a  SCTMS,  a  lumped  mass  representation  is  used,  i.e. 

H  dMp  =  Qp  dt  [20] 

where  Mp  =  Hass  of  PCM 

t  =  Time  of  storage 

H  =  Latent  heat  of  fusion 

Qp  =  Heat  load  of  PCM 

Solving  the  Eq.  20  for  the  differential  mass  of  the  PCM  yields 

dMp  =  0p  dt  /  H  [21] 

Since,  the  sensible  heat  storage  of  the  PCM  has  been  neglected,  due  to  the 
dominate  latent  heat  storage  of  the  PCM  the  above  equation  is  valid  for 
limited  cases  for  small  Stefan  numbers,  i.e. 

Ste  =  Cp  dTp  /  H  «  1  [22] 

where  Cp  =  Specific  heat  of  PCM 

dTp  =  Temperature  swing  of  PCM 
H  =  Heat  of  fusion  of  PCM 

In  addition,  the  thermal  resistance  to  conductive  heat  diffusion  in  the 
PCM  is  negligible  in  comparison  to  the  overall  coefficient  of  heat  transfer 
between  the  PCM  and  the  heat  transport  medium.  Therefore  Eq.  21  is  valid  for 
small  Biot  number,  i.e. 

Bi  =  UpARp  /  Kp  «  1  [23] 

where  Up  =  Overall  surface  coefficient  of  heat  transfer  in  the 
PCM  heat  exchanger 

A  Rp  =  Characteristic  thickness  of  the  PCM  layer 
Kp  =  Thermal  conductivity  of  the  PCM 

c.  Radiator:  To  determine  the  weight  of  the  radiator  in  a  SCTMS,  two 
principles  are  used,  a  mass  definition  and  a  heat  load  definition,  i.e. 


m 


dMr  =  pr  Lr  dAr 

[24] 

and 

dQr  =  qr  dAf 

[25] 

where  Mr  =  Mass  of  the  radiator 

Qr  =  Heat  load  of  the  radiator 
Ar  =  Radiative  surface  area 
pr  =  Mass  density  of  the  radiator 
Lr  =  Characteristic  thickness 
qf  =  Radiative  heat  flux 

Substituting  the  differential  area,  dAr,  from  Eq.  25  into  Eq.  24  yields 

dMr  =  (Pr  Lr  /  qr)  dQr  [26] 

Multiplying  and  dividing  Eq.  26  by  the  radiative  surface  area,  Ar  .g^s 
dM^  =  dQr  /  (  qr  Ar  /  pr  Lr  Ar) 

dMr  =  dQr  /  (  Qr  /  Mr) 

dMr  =  dQr  /  (  Gr)  [27] 

Hence,  if  the  specific  weight  of  the  radiator  can  be  assumed  to  be  a  constant, 
then  integrating  Eq.  27  yields  a  linear  relationship  between  radiator  mass  and 
thermal  load,  i.e. 

Mr  =  Qr  /  Gr  [28] 

B.  KNOWLEDGE-BASED  SYSTEM 

In  the  design  of  any  knowledge-based  system,  one  of  the  most  difficult 
problems  to  be  tackled  is  that  of  knowledge  representation  [14].  In  order  to 
first  understand  and  then  solve  a  problem,  one  must  be  able  to  represent  the 
objects,  relationships  and  the  abstract  concepts  of  the  domain  in  the 
artificially  intelligent  system.  Without  proper  representation  it  would  be 
difficult  to  achieve  a  significant  amount  of  intelligence  in  the  program,  as 
the  crux  of  'understanding'  a  problem  is  in  abstracting  it  to  a  representation 
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which  is  complete  and  amenable  to  manipulation  and  trans formation  that  leads 
to  a  solution.  Therefore,  the  initial  thrust  of  the  present  research  is 
directed  towards  the  development  of  appropriate  representation.  While 
addressing  the  representation  and  system  design  problems,  full  consideration 
is  given  to  the  state-of-the-art  in  the  field  of  knowledge-based  systems. 

1.  Automation  of  the  Design  Process:  A  typical  design  process  of  a  SCTMS 
includes  at  least  three  distinct  stages:  a  conceptual  design  stage,  a 
development  design  stage  and  a  design  verification  stage.  In  the  conceptual 
design  stage  system  tradeoffs  are  considered  prior  to  selecting  a  preliminary 
design.  The  thermal  design  of  SCTMS  requires  ‘systems  engineering'  since  a 
broad  spectrum  of  technical  disciplines  are  considered  in  this  design  process. 
As  spacecraft  systems  grow  in  size,  power  and  system  complexity,  the  number  of 
design  options  for  SCTMS  becomes  enormous.  Simple  tradeoff  studies  only 
consider  the  interdependence  of  the  various  design  parameters.  Many  current 
preliminary  design  processes  only  consider  a  single  point  design;  however, 
automation  of  the  preliminary  design  process  would  make  it  possible  to 
consider  total  payload  and  mission  requirements. 

To  assist  in  this  design  process,  a  knowledge-based  system  has  been 
developed  in  'Prolog'.  This  software  serves  as  a  design  assistant  for  the 
development  of  SCTMS.  The  design  assistant  is  implemented  according  to  the 
suggestions  found  in  [10].  These  suggestions  give  the  stages  of  operation  and 
the  levels  of  knowledge.  The  topmost  level  of  the  software  determines  an 
approximate  design  for  the  SCTMS.  The  lower  levels  then  in  turn  refine  the 
design,  within  their  specific  areas  of  expertise.  Fig.  10  presents  a  flow 
chart  of  the  overall  system  and  Fig.  11  illustrates  the  stages  of  operation 
and  the  levels  of  knowledge  required. 
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Specifically,  the  topmost  level  selects  a  PCM  which  is  compatible  with 
the  initial  temperature  specifications  and  determines  its  size  according  to 
the  mission  requirements.  This  process  is  facilitated  by  the  automatic  search 
and  selection  capabilities  of  PROLOG.  The  second  level  of  design  is  to  select 
a  heat-pipe,  radiator  combination  which  satisfies  the  given  mass  limit,  yet 
has  the  required  heat  transfer  capabilities.  This  combination  heat  pipe,  PCM 
and  radiator  then  serves  as  a  candidate  conceptual  design. 

Although  the  basic  heat  transfer  capabilities  of  the  selected  components 
are  matched,  the  actual  behavior  of  the  components  in  a  pulsed  power  system 
may  be  very  erratic.  Hence,  the  third  level  of  design  is  to  simulate  their 
behavior  under  the  mission  requirements,  using  the  transient  heat  flow  model 
described  above,  and  to  refine  the  overall  heat  transfer  coefficients  of  the 
components  so  that  the  system  behavior  pulse  after  pulse,  is  bounded  and 
predictable.  The  result  is  components:  heat  pipe,  PCM  and  radiator,  which 
together  satisfy  the  basic  mission  requirements. 

The  final  stage  of  the  design  is  to  determine  the  degree  of  enhancement 
(or  relaxation)  required  to  implement  components  with  those  heat  transfer 
performance  characteristics.  In  some  cases  of  PCMs,  heat  transfer  surface 
augmentation  is  necessary  to  compensate  for  relatively  low  thermal 
conductivity  of  the  selected  PCM.  In  other  cases  of  PCMs,  less  stringent 
requirements  of  the  overall  heat  transfer  coefficient  may  allow  for 
multi-configuration  options. 

2.  Knowledge-Base:  The  knowledge-base  of  the  design  assistant  has  three 
partitions.  Heat  Pipe,  PCMs,  and  Radiators.  To  build  this  knowledge-base  an 
extensive  literature  survey  covering  topics  of  latent  thermal  energy  storage 
systems  and  spacecraft  thermal  management  systems  was  made.  The  bibliography 
was  prepared  with  microcomputer  software  'The  Research  Assistant'  from  Hi-Q 
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(NAME )  (CLASS)  (TEMPERATURES  °c) 

(h*  -  KJ  /  KG>  (Q,.  Q,,  -  KJ  /  KG  /  C>  (k.  k»  - 
(DENSITIES  -  KG  /  M3)  (comment))) 


(polyethylene)  (organic  compound  polymer)  (132  132) 

(184)  (n-a  n-a)  (5. 18  4.6) 

(950950)  (solid-solid))) 

(polypropylene)  (organic  compound  polymer)  (166  166) 

(66.9)  (n-a  n-a)  (1.17  1.17) 

(901  901)  (solid-solid))) 

(polymethylene  oxide)  (organic  compound  polymer)  (181  181) 

(158.9)  (n-a  n-a)  (2.3  2.3) 

(1420  1420)  (solid-solid))) 

(polyethylene  oxide)  (organic  compound  polymer)  (57  57) 

(94.5)  (n-a  n-a)  (n-a  n-a) 

(1150  1150)  (sol  id-sol  id) > ) 

(pol y i soprene)  (organic  compound  polymer)  (32  35) 

(27.8)  (n-a  n-a)  (n-a  n-a) 

(n-a  n— a)  (solid-solid))) 

(85  7.  NaN03  -  8.8  7.  Na2S04  -  5.7  7.  NaCl  )  (inorganic  eutectic  salt)  (278  278) 
(168)  (n-a  n-a)  (0.78  0.78) 

(2290  2290)  ( i nexpensi ve ) > ) 

(95.4  7.  NaN03  -  4.6  7.  NaCl)  (inorganic  eutectic  salt)  (297  297) 

(178)  (n-a  n-a)  (0.69  0.69) 

<2260  2260)  (inexpensive))) 

(69.9  7.  CaC12  -  30.1  7.  NaCl)  (inorganic  eutectic  salt)  (490  490) 

(233)  (n-a  n-a)  (4  4) 

(2400  2400)  (inexpensive))) 

(46.2  7.  KC1  -  29.9  7.  CaCl  2  -23.9  7.  NaCl)  (inorganic  eutectic  salt)  (515  515) 
(265)  (n-a  n-a)  (3.8  3.8) 

(2160  2160)  (inexpensive))) 

(79  7.  LiOH  -21  7.  LiF)  (inorganic  eutectic  salt)  (427  427) 

(786)  (n-a  n-a)  (2.6  2.6) 

(1580  1580)  (high-storage) ) ) 

(LiOH)  (inorganic  compound  salt)  <460  460) 

(874)  (n-a  n-a)  (2.4  2.4) 

(1430  1430)  (hi gh-stor age) ) > 

(68  7.  KF  -  32  7.  LiF)  (inorganic  eutectic  salt)  (492  492) 

(523)  (n-a  n-a)  (3.5  3.5) 

(2520  2520)  (high-storage))) 


Figure  13.  PCM  Database 


(NAME)  (CLASS)  (TEMPERATURES  °c) 

(h«  -  KJ  /  KG)  <<^,_  Q,,  -  KJ  /  KG  /  C>  (k.  k»  -  W  /  M  /  C) 
(DENSITIES  -  KG  /  M}>  (coment))) 


(67.76  7.  A1  -  27  X  Cu  -5.25  X  Si)  (inorganic  eutectic  alloy)  (524  524) 
(308)  (n-a  n-a>  (256  256) 

(3300  3300)  ())) 

(66.8  7.  fll  -  33.2  7.  Cu)  (inorganic  eutectic  alloy)  (548  548) 

(275)  (n-a  n-a)  (273  273) 

(3510  3510)  () ) ) 

<82.1  7.  A1  -  13.3  7.  Si  -  4.6  7.  Mg)  (inorganic  eutectic  alloy)  (559  559) 
(423)  (n-a  n-a)  (199  199) 

(2590  2590)  ())) 

(87.4  7.  A1  -  12.6  X  Si)  (inorganic  eutectic  alloy)  (577  577) 

(431)  (n-a  n-a)  (180  180) 

(2661  2661)  <>>) 

(octadecane)  (organic  compound  other)  (28  28) 

(242.2)  (2.18  2.18)  (0.15  0.15) 

<910  770)  () ) ) 

(lithium  nitrate  trihydrate)  (inorganic  compound  hydrate)  (30  30) 

(296.5)  (2.09  2.09)  (0.5  0.5) 

(1550  1430)  () ) ) 

(Na2HPQ4  «  12  (H20) )  (inorganic  compound  hydrate)  <36  36) 

(279.8)  (1.67  1.93)  (0.51  0.48) 

<1520  1450)  ())) 

(myristic  acid)  (organic  compound  acid)  (58  58) 

(200.5)  (1.59  2.26)  (0.1  0.1) 

(860  860)  () ) ) 

(cerrobend)  (inorganic  compound  alloy)  (70  70) 

(33.4)  (0.17  0.17)  <18.8  18.8) 

(9400  9200)  () ) ) 

(acetimide)  (organic  compound  other)  (81  81) 

(242.2)  (2.8  2.8)  (0.25  0.25) 

(1160  n-a)  () ) ) 

(CaC12  *  6  (H20) )  (inorganic  compound  hydrate)  (30  30) 

(170.3)  (n-a  n-a)  (108.73  50.97) 

(1802  1562)  (incongruent) ) ) 

(Zn  (N03)  2*6  (H20) )  (inorganic  compound  hydrate)  (36.1  36.1) 

(133.9)  (n-a  n-a)  (n-a  46.83) 

(1937  1828)  (expensive))) 


Figure  13.  continued 


W  /  M  /  C> 


(NAME)  (CLASS)  (TEMPERATURES  °c> 

(h,  -  KJ  /  KG)  (Q,_  -  KJ  /  KG  /  C)  (k.  k,  - 

(DENSITIES  -  KG  /  M3)  (comment))) 


(Ca  (NQ3)  2*4  (H20) )  (inorganic  compound  hydrate)  (47.2  47.2) 

(161.9)  (n-a  n-a>  (n-a  n-a) 

(n-a  n-a)  < ) > ) 

(Mg  (N03)  2*6  (H20) )  (inorganic  compound  hydrate)  (90  90) 
(159)  (n-a  n-a)  (61.5  48.95) 

(1636  1550)  (corrosi ve) ) ) 

(MgC12  *  6  <H20) >  (inorganic  compound  hydrate)  (115  115) 

(164.9)  (n-a  n-a)  (69.45  56.87) 

(1569  1450)  (corrosive))) 

(wax  suntech  pll6)  (organic  eutectic  paraffin)  (44  44) 

(226)  (2.95  2.51)  (0.24  0.24) 

(818  760)  (nonmonotonic-conductivity) ) ) 


Figure  13.  continued 


Table  V.  Parameters  Considered  in  the  Analysis 

*  Pressure  drop  through  packed  bed 

*  Flow  rate  of  HTF 

*  Test  section  geometry 

*  Total  heat  stored 

*  Melt  time 

*  Effect  of  cycling 

*  Convective  heat  transfer  coefficient  between  HTF  and  PCM 

*  Porosity 
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Microware,  #8-13322 ,102A  Avenue,  Surrey,  B.C.  Canada  V3T  5J7.  The  literature 
is  organized  under  three  headings,  SCTMS,  latent  thermal  energy  storage  and 
phase  change  mathematical  models. 

a.  Phase  Change  Materials  Knowledge-Base:  This  consists  of  properties  of 
about  50  selected  PCMs.  The  PCMs  are  grouped  into  different  classes  as 
shown  in  the  Fig.  12.  The  PCMs  were  selected  keeping  in  view  the  desired 
properties  shown  in  the  Table  III  and  IV.  A  sample  PCM  database  is  shown  in 
the  Fig.  13.  The  columns  contain  the  name  of  the  PCM,  melting  point,  heat  of 
fusion,  thermal  conductivity  in  solid  state  and  liquid  state,  specific  heat  in 
solid  state  and  liquid  state  and  the  last  column  contains  comments  with 
respect  to  PCMs  considerations  shown  in  the  Table  V. 

b.  Heat  Pipe  Knowledge-Base:  At  present,  the  heat  pipe  knowledge-base 
consists  of  single  point  database.  The  heat  pipe  is  chosen  according  to  the 
specific  weight.  The  creation  of  a  more  complete  heat  pipe  database  will 
involve  a  greater  effort  due  to  its  complex  nature  and  all  the  factors 
involved  are  not  known  at  this  stage.  Figure  14  gives  some  ideas  in  that 
direction.  The  implementation  of  a  more  complete  knowledge-base  for  heat 
pipes  may  be  performed  as  a  future  task. 

c.  Radiator  Knowledge-Base:  Here  again,  a  single  database  is  considered 
for  the  radiator  knowledge-base.  The  radiator  is  selected  according  to  the 
specific  weight.  Again,  similar  to  the  heat  pipes,  the  creation  of  a  more 
complete  database  for  the  radiator  has  been  deferred  to  a  future  task. 


IV.  EXPERIMENTAL  PROGRAM 

As  an  ongoing  research  task,  this  report  primarily  addresses  properties 
and  thermal  behavior  of  selected  PCMs  as  candidates  for  possible  use  as 
thermal  buffers  of  burst  power  conditions  of  SCTMS  materials. 


A.  FUSIBLE  PELLET  CONCEPT 


In  the  present  concept,  a  packed  bed  of  pellets  can  be  designed  to 
provide  a  specified  thermal  energy  storage  and  charging/discharging  rates.  The 
two  PCMs  (Calcium  choloride  hexahydrate  and  high  density  polyethylene) 
selected  for  evaluation  have  the  potential  to  satisfy  the  requirement  of  being 
impervious  to  water,  the  heat  transfer  fluid  used  in  the  experimental  tests. 

In  the  case  of  the  calcium  chloride  hexahydrate,  the  pellets  or  capsules  are 
small,  self-contained  units  in  which  the  active  core  is  separated  from  the 
water  by  organic  resins  making  it  impervious  to  the  water.  While  in  the  case 
of  the  High  Density  PolyEthylene  (HOPE),  the  PCM  is  non-porous  and  no  coating 
or  encapsulation  is  necessary. 

The  design  and  analysis  of  packed-bed  heat  exchangers  require  both 
thermophysical  and  geometric  information,  such  as  the  porosity,  device 
length-to-diameter  ratio,  and  numerous  thermodynamic  properties  of  the  thermal 
storage  material.  In  a  preliminary  research  task,  some  of  these  properties 
were  obtained.  Special  attention  was  given  to  the  HDPE  because  of  its 
apparent  unique  characteristics.  In  the  recent  work  on  solid-solid  PCMs,  most 
of  the  effort  has  been  directed  at  the  solid-solid  transition  and  how  to 
produce  the  PCM,  such  as  silanegrafted  crosslinking  or  electron  beam 
irradiation.  Relatively  little  work  has  been  devoted  to  the  other  properties 
of  the  material.  However,  few  research  teams  have  addressed  the  problems  of 
the  HDPE,  such  as  its  compressibility,  relatively  low  thermal  conductivity  and 
compatibility  with  different  heat  transfer  fluids  [16-18]. 

B.  HEAT  TRANSFER  EXPERIMENTS 

The  development  of  packed  bed  heat  exchangers  for  SCTMS  requires  new 
technology.  Emphasis  is  placed  on  developing  a  method  for  predicting  the 


thermal  response  of  an  entire  packed  bed  heat  exchanger  unit  rather  than  a 
single  PCM  pellet,  in  order  to  facilitate  its  design. 

All  of  the  existing  latent  heat  storage  devices  use  sensible  heat 
transfer  and  suffer  significant  losses  due  to  the  Second  Law  of 
Thermodynamics.  These  losses  can  be  minimized  if  a  latent  heat  transfer 
process  is  used,  i.e.  a  boiling  and  condensing  heat  transfer  fluid.  Hence  the 
idea  of  a  packed  bed  with  a  PCM  integrated  into  the  condenser  in  a  heat  pipe 
is  promising.  However,  if  one  selects  steam  as  the  heat  transfer  fluid, 
several  problems  with  PCMs  have  been  observed  such  as  pellet  adhesion  and 
moisture  adsorption.  For  the  crosslinked  HDPE  with  a  melting  point  of  about 
130  °C,  ethylene  glycol,  silicone  oil,  water  and  steam  have  been  tested.  To 
test  the  process,  the  heat  exchanger  in  Figure  15  was  used.  It  consists  of  a 
rectangular  cross  section,  stainless  steel  column  to  hold  the  pellets,  and  has 
a  rectangular  glass  window  on  one  side  for  photographic  observation  of  the 
phase  change  phenomenon  (Fig.  16).  The  design  and  analysis  of  packed  bed  heat 
exchanger  requires  both  thermophysical  and  geometric  information.  Table  V. 
shows  the  important  parameters  considered  in  the  analysis.  Both  pressure  drop 
and  temperature  measurements  of  the  heat  transfer  fluid  are  being  considered. 

C.  DESIGN  ANALYSIS 

In  order  to  analyze  the  problem  of  using  PCM  pellets  in  a  packed  bed  heat 
exchanger,  several  simplifying  assumptions  have  been  made.  They  are  as 
follows: 

a.  Heat  transfer  fluid  properties  are  assumed  constant, 

b.  Heat  transfer  coefficient  between  the  heat  transfer  fluid  and 
the  PCM  is  uniform  in  the  heat  exchanger, 

c.  The  mean  velocity  of  the  heat  transfer  fluid  in  the  axial 
direction  is  constant. 


d.  The  axial  conduction  heat  transfer  between  the  pellets  is 
negligible  compared  with  the  heat  transfer  via  the  fluid, 

e.  The  temperature  distribution  of  the  PCM  is  one-dimensional  in 
the  axial  direction, 

f.  The  pellets  are  uniform  in  size  and  spherical  in  shape, 

g.  The  inlet  temperature  of  the  heat  transfer  fluid  is  constant. 

The  conservation  of  energy  in  the  packed  bed  heat  exchanger  can  be 

described  by  the  following  relationship  between  the  time  rate  of  change  of  the 
PCM  and  the  convective  coefficient  of  heat  transfer: 

hcAT  As  =  dMs/dt  H  [29] 

where  H  =  Enthalpy  of  the  PCM. 

Integrating  with  respect  to  time  from  the  initial  time  of  phase  change  to 
the  time  when  one  pellet  has  phase  changed  yields: 

x  =  Ms  H  /  (hcAT  As)  [30] 

To  size  the  length  of  a  packed  bed  heat  exchanger,  one  can  check  the  length 
that  the  fluid  flows  during  the  time  for  one  pellet  layer  to  melt.  An 
expression  can  be  obtained  for  this  length  as  follows: 

<5  =  4  Ps  H  Ds  V  /  (6  hcATTre  D  )  [31] 

Following  the  analysis  of  T.  Asahina  and  M.  Kosaka  [19],  the  pressure  drop  of 
the  fluid  passing  through  the  packed  bed  can  be  obtained  from  the  following 
expression: 

AP  =  (2  fpL  V2  /  (g  Ds  Acs) )  (1  -e)  /  e  [32] 

where  f  =  1.75  /  Re®** 

and  Re  =p0s  /  u  (4V  /  (D2  (1  -  e  2/ 3) ) 
for  Re  >  200 

Table  VI  gives  several  properties  of  the  two  PCMs  pellets  being 
experimentally  investigated.  The  calcium  chloride  hexahydrate  (CaCl2.6H20) 
material  is  being  tested  first  in  a  rectangular  cross  section  packed  bed  since 


Table  VI.  Properties  of  PCM  Pellets 
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this  arrangement  allows  for  easy  photographic  recording  of  the  phase  change 
phenomena. 

V.  EXPERIMENTAL  RESULTS 

From  an  array  of  potential  materials,  two  candidate  PCMs  in  pellet  form 
were  selected  for  review:  encapsulated  calcium  chloride  hexahydrate  and 
radiation-crosslinked  high-density  polyethylene  (HDPE).  The  radiation 
cross-linking  gives  the  HOPE  the  unique  physical  property  of  being  self 
encapsulated.  The  HDPE  can  be  heated  to  temperatures  in  excess  of  150  °C 
without  becoming  entirely  liquid.  The  individual  pellets  of  HDPE  are 
form-stable.  Recent  studies  [15]  have  shown  that  the  HDPE  pellets  of  retain 
their  shape;  however,  pellets  do  tend  to  stick  together.  Lane's  recent  book 
[16]  on  latent  heat  storage  materials  gives  an  excellent  and  comprehensive 
review  of  both  of  these  materials. 

In  order  to  study  the  direct  contact  heat  transfer  between  the 
form-stable  HDPE  or  the  encapsulated  CaCl2.6H20  and  the  heat  transfer  fluid, 
several  material  properties  were  needed  for  the  design  of  the  packed  bed  heat 
exchanger  including  the  porosity,  the  pellet  microhardness  (Vickers  Hardness 
Number),  the  melting  temperature  and  the  recrystallization  temperature. 

A  recent  investigation  [17]  has  shown  that  HDPE  has  a  melting  point  in 
the  temperature  range  of  125-140  °C  and  values  for  latent  heat  as  high  as  210 
kj/kg  which  we  have  confirmed  experimentally.  HDPE  melts  and  freezes 
reversibly  and  congruently.  Thus  the  value  of  heat  of  fusion  remains  at  a 
constant  value  even  after  repeated  cycling  36  through  the  melting  point.  HDPE 
changes  from  white  opaque  crystalline  solid  to  a  clear  transparent 
viscoelastic  liquid  phase.  A  finely  divided  white  residue  was  observed  to  be 
present  at  the  bottom  of  the  flask  when  heated  to  a  temperature  higher  than 
the  melting  point  in  a  bath  of  water  or  oil.  It  has  been  observed  [18]  that 
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HDPE  pellets  which  are  exposed  to  a  dose  of  about  8  megarads  adhere  to  one 
another  at  their  point  of  contact,  but  the  packed  bed  of  pellets  maintain  its 
open  porous  structure  throughout  the  thermal  cycling  phase. 

1.  Porosity  and  Density:  The  porosity  of  the  packed  bed  heat  exchanger 
was  calculated  from  experimental  measurements  of  the  HDPE  and  calcium  chloride 
hexahydrates.  The  pellets  were  poured  into  a  calibrated  glass  cylinder  while 
varying  packing  pressure  was  axially  applied  with  a  dead-weight  system.  The 
resulting  packing  pressure  was  measured  on  a  calibrated  load  cell  with  an 
analog  readout.  Water  was  added  to  the  packed  bed  and  measured  to  determine 
the  void  volume  as  a  function  of  both  bed  temperature  and  packing  pressure. 

The  temperature  of  the  packed  bed  in  the  cylinder  was  maintained  by  a  small 
tape  heater.  In  addition,  a  microhardness  tester  was  used  to  determine  the 
hardness  of  the  HDPE  as  a  function  of  radial  location  in  the  pellet.  During 
the  cross-sectioning  of  an  encapsulated  bed  of  HDPE  for  the  microhardness 
test,  a  small  void  in  the  center  core  of  each  pellet  was  discovered.  Hence, 
the  mass  density  of  the  hollow  pellets  was  also  calculated  from  experimental 
measurements  of  the  HDPE  to  be  914.kg/m^.  To  eliminate  any  initializing  effect 
of  the  applied  packing  pressure,  multiple  loading/unloading  cycles  were  made 
until  consistent  values  of  porosity  were  obtained.  Table  VII  gives  the 
calculated  values  of  the  porosity  for  the  HDPE  as  a  function  of  packing 
pressure  and  bed  temperature.  Table  VIII  gives  the  calculated  values  of  the 
porosity  for  calcium  chloride  hexahydrate  as  a  function  of  packing  pressure  at 
the  melting  temperature. 


Table  VII.  HDPE  Porosity  versus  Temperature  and  Pressure 


Packing  Pressure  [kPa]  Bed  Temperature  in  degrees  centigrade 


60 

80 

100 

120 

103.4 

0.3386 

0.3400 

0.3428 

0.3370 

155.1 

0.3371 

0.3386 

0.3410 

0.3342 

106.8 

0.3355 

0.3368 

0.3382 

0.3329 

258.5 

0.3340 

0.3350 

0.3361 

0.3312 

Table  VIII.  CaCl2.6H20  Porosity  versus  Pressure 


Packing  Pressure  [kPa] 


0 

103.4 

258.5 


Bed  Temperature  =  27  degrees  centigrade 
Porosity 
0.4822 
0.4800 
0.4772 


2.  Photographic  Observations:  An  array  of  nine  HDPE  pellets  was  attached 
to  an  Aluminum  block  (0.057  m  x  0.057  m  x  0.006  m)  having  a  fine  grid  (0.0013 
m)  machined  on  the  surface.  Two  400  watt  cartridge  heaters  were  placed  in  two 
holes  drilled  into  the  side  of  the  aluminum  block.  Photographs  of  the  HDPE 

pellets  were  taken  at  temperatures  ranging  from  117  to  175  °C.  From  negatives 

(0.10  m  x  0.13  m)  of  the  photographs,  enlarged  prints  (0.20  m  x  0.26  m)  of  the 
pellets  were  made  to  determine  the  variation  of  size  and  shape  with  respect  to 

temperature.  The  resulting  prints  also  showed  that  during  the  phase  change 

process  the  HDPE  pellets  changed  from  opaque  to  transparent.  The  nominal 


diameter  of  the  pellets  was  measured  to  be  approximately  0.004  m.  Figure  16 

% 

shows  three  different  temperature  conditions  of  a  single  pellet.  The  right 
photograph  has  a  circle  superimposed  around  the  pellet  illustrating  the 
variation  from  an  idealized  spherical  shape.  The  middle  photograph  shows  the 
same  pellet  during  the  phase  change  with  an  opaque  center  and  a  transparent 
outer  shell.  The  left  photograph  in  Fig.  16  illustrates  a  transparent  pellet 
after  complete  phase  change  has  occurred.  The  photographs  were  taken  with  a 
Calumet  4x5  with  double  bellows  extension.  The  lens  was  a  F4. 8/0. 216  m 
Ilex.  The  shutter  speed  range  was  1  to  1/150  sec  and  F-stop  range  was  4.8  to 
45.  The  lighting  was  accomplished  with  two  3800  K  photo  spot  lights.  The 
photographs  were  taken  at  1/8  sec  and  F16  with  a  resulting  1.2  to  1  ratio  of 
the  negative  to  original  length.  The  film  selected  was  an  ASA  50  Polaroid 
positive/negative  4x5  land  film. 

Several  other  characteristics  of  the  HDPE  pellets  in  their  so-called 
'form-stable'  liquid  phase  were  observed.  A  pellet  at  175°C  was  squeezed  by  a 
pair  of  hand  tweezers.  The  pellet  was  noted  to  be  soft  and  pliable  and  showed 
a  strong  tendency  to  adhere  to  the  metal  tweezers.  This  adhesion  property  was 
reconfirmed  in  Fig.  17  which  shows  a  packed  bed  of  pellets  having  undergone 
cyclic  heating  and  packing  in  a  graduated  cylinder  during  a  porosity  test. 

The  pellets  clearly  adhere  to  each  other  even  when  heated  in  a  bath  of  oil. 

3.  HDPE  Pellet  Microhardness :  A  Leitz  MINILOAD-Hardness  Test,  equipped 
with  a  square  Vickers  pyramid  indenter,  was  used  to  measure  the  hardness  of 
the  crosslinked  HDPE  at  various  radial  locations  on  numerous  sectioned  pellets 
encapsulated  in  a  hard  resin.  The  Vickers  hardness  number  is  defined  as  the 
load  divided  by  the  area  of  contact  from  the  indenter.  The  results  of  16 
different  measurements  showed  that  the  Vickers  hardness  number  varied  from  5.8 
to  9.5  kg/mm2.  These  results  indicate  that  the  crosslinked  HDPE  is  a  very  soft 


material  in  comparison  to  metals.  The  general  trend  of  the  tests  indicated 
that  the  material  is  harder  near  the  center  core  and  becomes  softer  near  the 
outer  surface. 

The  microhardness  testing  was  restricted  to  the  HDPE  pellets  since  the 
salt  hydrates  undergo  a  normal  phase  change  process  with  the  liquid  being 
encapsulated  in  a  flexible  spherical  shell. 

4.  HDPE  Pellet  Differential  Thermal  Analysis  Results;  To  determine  the 
melting  and  the  recrystal lization  temperatures ,  a  thermogal vanometry- 
differential  thermal  analysis  (TG-OTA)  was  performed.  The  TG-DTA  was 
originally  the  only  technique  available.  Its  capability  is  limtied  to 
detection  of  phase  transitions  and  chemical  reactions.  The  TG-OTA  is  not  to 
be  confused  with  the  Thermogravimetric  analysis  (TGA)  used  for  composition  and 
thermal  stability  studies.  Several  random  samples  of  the  pellets  were  used. 

The  results  indicate  that  the  melting  temperature  of  the  crosslinked  HDPE  is 
approximately  130  °C  and  the  recrystallization  temperature  is  approximately 
120  °C.  Differential  scanning  calorimetry  (DSC)  results  would  have  allowed 
specific  heats  of  the  crosslinked  HDPE  to  be  measured.  However  due  to 
nonavailability  of  a  DSC  at  this  facilities,  the  specific  heats  were  not 
determined.  Others  [17]  have  recently  reported  such  results  for  crosslinked 
HDPE.  It  should  be  noted  that  the  values  from  the  DSC  are  dependent  on  the 
thermal  history  due  to  the  nature  of  polymers.  Hence,  one  might  conclude  that 
actual  test  of  prototype  packed  bed  heat  exchangers  are  necessary.  Hence, 
actual  testing  of  a  prototype  packed  bed  heat  exchanger  is  planned  for  several 
PCH's. 

5.  Parametric  Studies:  Figures  18-21  give  the  results  of  a  parametric 
study  of  the  pressure  drop  as  a  function  of  the  ratio  of  the  length  to 


diameter  of  the  packed  bed.  A  set  of  lines  for  constant  diameters  are  drawn 
for  a  fixed  volumetric  flow  rate  of  the  heat  transfer  fluid  resulting  in 
constant  Reynolds  number  lines.  Water  was  selected  as  the  heat  transfer  fluid 
in  the  calculations.  In  the  case  of  cylindrical  packed  beds  of  HDPE,  one  can 
see  in  Fig.  18  that  for  a  Reynolds  number  equal  to  289  and  a  length  to 
diameter  ratio  of  6  that  the  pressure  drop  would  be  approximately  26  kPa. 
Likewise  for  a  CaCl^^O  filled  cylindrical  packed  bed,  one  can  see  in  Fig. 

19  that  for  a  Reynolds  number  equal  to  294  and  a  length  to  diameter  ratio  of  6 
that  the  pressure  drop  would  be  approximately  2.3  kPa.  In  the  case  of  a 
rectangular  cross  section  packed  bed  heat  exchanger.  Fig.  20  and  21  illustrate 
the  relationship  between  the  pressure  drop  and  the  length  to  hydraulic 
diameter  ratio  for  HDPE  and  CaCl2.6H20,  respectively.  It  is  evident  from  Eq. 
31 

6  =  4  ps  H  D$  V  /  (6  hc  ATttc  Dp), 

that  the  pellet  diameter  along  with  the  porosity  influence  the  pressure  drop 
across  the  packed  bed  heat  exchanger. 

6.  Thermal  Analysis:  Thermal  analysis  of  a  packed  bed  of  CaCl2.6H20  in 
pellet  form  was  performed.  Water  was  chosen  for  the  Heat  Transfer  Fluid  (HTF) 
as  it  is  clear,  abundantly  available,  has  suitable  properties  and  is 
inexpensive.  For  this  experimental  application  (see  Fig.  22),  the 
availability  of  low  pressure  steam  and  tap  water  enabled  a  variety  of  flow  and 
temperature  conditions  to  be  obtained  by  proper  mixing.  The  storage  unit 
instrumentation  consisted  of  a  paddle  wheel  flowmeter,  weight  tank,  a 
thermometer,  two  pressure  gages  and  eight  thermocouples.  The  paddle  wheel 
flowmeter  and  weight  tank  was  used  to  measure  the  water  flow  rate  and 
thermocouples  measured  temperatures  at  various  locations  in  the  storage  unit. 
These  were  used  to  determine  whether  the  PCM  was  charged  or  discharged  at  the 
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beginning  of  each  test  and  to  obtain  a  qualitative  understanding  of  the 
thawing  or  freezing  process. 

Tests  were  conducted  for  both  charging  and  discharging  of  the  PCM.  Prior 
to  starting  each  test,  the  PCM  was  cooled  or  heated  so  that  it  was  completely 
frozen  for  charging  tests  and  melted  for  discharging  tests.  The  temperature 
of  the  PCM  was  verified  by  reading  the  thermocouples  to  ensure  that  all  PCM 
temperature  readings  were  at  least  3°F  below  the  melting  point  for  charging 
process  and  3°F  above  the  melting  point  for  discharging  process. 

Six  test  runs  were  performed,  three  charging  and  three  discharging,  with 
water  flow  rates  of  0.2,  0.4  and  0.6  gallon  per  minute.  The  test  conditions 
were  established  by  operating  the  steam  mixing  valves  and  the  valve  on  the 
water  line.  The  data  acquisition  for  discharging  process  began  once  all  the 
thermocouples  readings  were  approximately  84°F.  The  inlet  HTF  temperature  was 
lowered  to  approximately  76°F  by  shutting  off  the  steam  and  the  temperatures 
and  flow  rates  were  recorded  at  an  interval  of  30  seconds  until  the 
discharging  process  was  completed.  Figures  23,  24  and  25  show  the  temperature 
vs.  time  curves  for  the  mass  flow  rates  0.2,  0.3  and  0.4  gpm.  For  the  three 
mass  flow  rates  the  outlet  temperature  lags  behind  in  inlet  temperature, 
showing  that  the  PCM  is  discharging  heat.  It  is  difficult  to  measure  the  time 
required  to  melt  the  PCM,  as  the  heat  capacity  of  the  heat  exchanger  is  far 
greater  than  the  latent  heat  of  PCM.  Hence  the  response  of  the  PCM  packed  bed 
cannot  be  accurately  determined.  A  steady  state  analysis  cannot  be  performed 
since  it  is  difficult  to  lower  the  inlet  water  temperature  from  86°F  to  76°F 
instantaneously  with  the  present  setup.  A  transient  analysis  was  performed  to 
approximately  calculate  the  mass  of  PCM  discharged  vs.  time  for  different  mass 
flow  rates  using  the  following  equation. 

Considering  an  energy  balance  at  any  instant  of  time,  the  heat  absorbed 
by  the  water  plus  the  heat  lost  to  the  surrounding  due  to  natural  convection 
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Figure 225.  Present  Experimental  Setup 
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equals  the  heat  discharged  by  the  PCM, 
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where  m 


water 


=  mp  X  (t  )  H 

=  mass  flow  rate  of  water 


mp  =  mass  of  PCM 

X  =  fraction  of  PCM  melted 

Tout  =  ou^et  water  temperature 

Tin  =  inlet  water  temperature 

h  =  film  coefficient  of  water 

A  =  surface  area  of  the  heat  exchanger 

Tpox  =  temperature  of  the  surface  of  the  heat  exchanger, 

'•e-  lTin  *  Tout>  '  2-° 

T$ur  =  ambient  temperature 

t  =  time 

The  sensible  heat  contribution  is  small  compared  with  the  latent  heat  of 
the  PCM.  Figure  26  shows  mass  fraction  vs  time  curves  for  different  flow 
rates.  The  time  required  for  discharging  the  PCM  increases  as  the  flow  rate 
drops.  These  results  should  be  used  only  as  a  guide  to  future  experimental 
work.  Some  recommendations  for  future  experimental  work  are  discussed  in  the 
next  section. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  application  of  phase  change  materials  for  thermal  buffering  of 
temperature  sensitive  equipment  from  pulsed  heat  loads  of  future  SCTMS  appears 
to  satisfy  the  objectives.  The  initial  selection  of  the  candidate  materials 
might  be  ranked  on  latent  heat  of  fusion,  however,  the  final  selection  should 
will  not  be  based  solely  on  a  maximized  latent  heat  of  fusion.  The  specific 
weights  of  radiator  and  the  heat  pipes  are  operationally  dependent  on  numerous 
parameters.  Novel  heat  pipe  and  radiator  designs  will  be  required  to  achieve 
the  desired  performance. 

The  results  of  this  preliminary  investigation  of  the  PCMs  for  SCTMS  have 
been  presented  to  identify  the  feasibility  and  capabilities  of  the  selected 
PCMs  for  thermal  buffers  in  fusible  pellet  concept.  Several  technical  issues 
relating  to  the  thermophysical  properties  were  evaluated  to  determine  the 
design  improvements  for  the  experimental  facilities.  The  highly  crystalline 
polymers  such  as  radiation  crosslinked,  high  density  polyethylene  show 
potential  for  a  direct  heat  transfer  process  between  a  boiling/condensing  heat 
transfer  fluid  and  a  PCM  due  to  its  rather  unique  property  of  self 
encapsulation.  Further  experimental  verification  of  the  bulk  thermophysical 
properties  of  pelletized  PCMs  in  the  packed  bed  form  is  needed.  Also  methods 
to  avoid  adhesion  of  HOPE  pellets  at  temperature  up  to  150°C  for  packed  bed 
conditions  are  needed. 

The  results  of  the  thermal  analysis  of  the  PCM  storage  unit  recommends  an 
altogether  new  experimental  set  up.  The  proposed  experimental  setup  is  shown 
in  Fig.  27.  It  consists  of  the  following: 


a. 

High  temperature  bath 

b. 

Bypass  valves 

c. 

Air  cooler 

d. 

Thermocouples 

e. 

PCM  heat  exchanger 

f. 

Micro-computer  data 

acquisition  system 

By  having  this  arrangement  the  following  problems  associated  with  the 
current  setup  are  resolved. 

a.  Control  of  the  input  fluid  temperature. 

b.  Sensible  heat  of  the  heat  exchanger  minimized. 

c.  Fluid  leakage. 

d.  Lower  fluid  pressures  in  the  system  by  using  oil  as  HTF  for 
higher  temperature  approximately  instead  of  steam. 

The  prototype  design  assistant  satisfies  our  initial  expectations.  It 
properly  selects,  sizes,  and  refines  the  components  for  a  three  component 
SCTMS  model.  Our  next  goal  is  to  install  a  larger  database  of  components  and 
mission  requirements  and  to  exercise  the  software  over  a  larger  set  of  design 
problems. 

Two  interesting  extensions  of  this  work  are  possible.  The  first  is  to 
extend  this  work  horizontally:  to  develop  a  family  of  design  assistants,  each 
tailored  to  a  specific  type  of  design  problem,  continuous  power,  very  high 
power  pulsed  systems,  pulsed  systems  with  very  short  cycle  time,  etc.  These 
separate  assistants  could  then  be  called  by  a  supervisory  process  which  would 
determine  from  the  mission  requirements  which  assistant  should  be  consulted. 
The  second  extension  would  be  vertical:  to  develop  additional  levels  to  the 
present  assistant  which  would  further  refine  the  physical  dimensions  and 
structural  characteristics  of  a  three  component  SCTMS. 

A  very  interesting  prospect,  from  an  AI  perspective,  is  the  possibility 
of  capturing  the  knowledge  of  SCTMS  design,  as  it  is  acquired  by  the 
researchers  themselves,  documenting  it,  implementing  it,  and  preserving  it  in 
a  body  of  software.  Hence  with  a  working  prototype  in  hand,  the  most 
important  phase  of  AI  work  begins.  A  significant  amount  of  domain  specific 
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